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CO 

o 
o 



A. Yu. Aladyshkin, W A. I. Buzdin, ^ A. A. Fraerman, W A. S. Mel'nikov, W D. A. Ryzhov, W* A. V. Sokolov W 

^' Institute for Physics of Micro structures, Russian Academy of Sciences, 603950, Nizhny Novgorod, GSP-105, Russia, 
^ Centre de Physique Moleculaire Optique et Hertzienne, Universite Bordeaux 1-UMR 5798, CNRS, F-33405 Talence Cedex, 

France 
(February 2, 2008) 

On the basis of phenomenological Ginzburg-Landau approach we investigate the problem of order 
^■y parameter nucleation in hybrid superconductor/ferromagnetic (S/F) systems with a domain struc- 

ture in applied external magnetic field. Both the isolated domain boundaries and periodic domain 
structures in ferromagnetic layers are considered. We study the interplay between the superconduc- 
tivity localized at the domain walls and far from the walls and show that such interplay determines 
a peculiar field dependence of the critical temperature T c . For a periodic domain structure the 
cn : behavior of the upper critical field of superconductivity nucleation near T c is strongly influenced by 

the overlapping of the superconducting nuclei localized over different domains. 
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I. INTRODUCTION 



The problem of coexistence of superconducting and magnetic orderings has been studied for several decades (see, 
e.g., Refs. 1 ' 2 for review). One can separate two basic mechanisms responsible for interaction of superconducting 
order parameter with magnetic moments in the ferromagnetic state: (i) the electromagnetic mechanism (interaction 
i-^j | of Cooper pairs with magnetic field induced by magnetic moments) which was first discussed by V. L. Ginzburg 3 
in 1956; (ii) the exchange interaction of magnetic moments with electrons in Cooper pairs. The revival of interest 
to the fundamental questions of magnetism and superconductivity coexistence has been stimulated, in particular, 
by the recent investigations of the hybrid superconductor/ferromagnetic (S/F) systems. Such thin film structures 
consist of a ferromagnetic insulator film and superconducting film deposited on it. Similar situation can be obtained 
with a metallic ferromagnet when a superconducting film is evaporated on the buffer oxide layer in order to avoid 
proximity effect. The superconducting properties of such structures attract a growing interest due to a large potential 
for applications. In particular, such hybrid S/F systems are intensively investigated in connection with the problem of 
controlled flux pinning. The enhancement of the depinning critical current density j c has been observed experimentally 
' for superconducting films with arrays of submicron magnetic dots, 4-6 antidots, 7 and for S/F bilayers with domain 
structure in ferromagnetic films. 8 Theory of vortex structures and pinning in the S/F systems at rather low magnetic 
^f*) ' fields (in the London approximation) has been developed in Refs. 9-18 . 

A nonhomogeneous magnetic field distribution induced by the domain structure in a ferromagnetic layer influences 
strongly the conditions of the superconducting order parameter nucleation, and, as a consequence, the hybrid S/F 
systems reveal a nontrivial phase diagram in an external applied magnetic field H (see, e.g., Refs. 19-21 ). In this 
paper we focus on the theoretical study of this phase diagram on the basis of phenomenological Ginzburg-Landau 
(GL) model. We assume that the electromagnetic mechanism mentioned above plays a dominant role and neglect 
the exchange interaction which is obviously suppressed provided superconducting and ferromagnetic layers are well 
separated by an insulating barrier. We also assume that the domain walls are well pinned and do not take account of 
£j , changes in the domain structure with an increase in H. 

The distribution of the magnetic field induced by the domain structure is determined by the ratio of two length 
scales: thickness of ferromagnetic film D and distance between the domain walls w (hereafter we neglect a finite width 
of the domain wall, i.e. consider this width to be much less than the superconducting coherence length). Provided the 
ferromagnetic film is rather thick (D 3> w), the magnetic field in a thin superconducting film is almost homogeneous 
over the domain and suppresses the critical temperature of superconductivity nucleation. In this case with the decrease 
in the temperature the superconductivity must firstly appear just above the domain wall (see Refs. 22,23 ) due to the 
mechanism analogous to the one responsible for the surface superconductivity below H c ^ (see Ref. 24 ). Thus, in this 
limit the domain walls stimulate the nucleation of the superconducting order parameter. Note, that the same effect 
should reveal for two-dimensional magnetic field distributions induced, e.g., by magnetic dots, and results in the 
dependence of the upper critical field on the angular momentum of the superconducting nucleus wave function (see 
Refs. 25-27 ). 

For a thin ferromagnetic film (D -C w) the magnetic field decays with the increase in the distance from the domain 
wall and almost vanishes inside the domain. In the absence of the external field such domain wall should locally 
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weaken superconductivity as it was discussed in Ref. 9 . The superconducting nucleus in this case should appear far 
from the domain wall. As we switch on an external magnetic field, we can control the position of the superconducting 
nucleus suppressing the order parameter inside the domains. Thus, the phase diagram of the S/F bilayer is generally 
determined by the interplay between the superconductivity nucleated at the domain walls and in between these walls. 
For small period domain structures (when w is comparable with the nucleus size) this simple physical picture based 
on consideration of isolated superconducting nuclei should be modified taking account of the interaction between the 
superconducting nuclei localized above different domain walls. 

Our further consideration is based on the linearized GL equation for the order parameter ^: 

-(v + £a)'.-^.. 

Here A(r) is the vector potential, B(r) = V x A(r), $o is the flux quantum, £(T) = £0/ \A ~~ T/T c q is the coherence 
length, and T c0 is the critical temperature of the bulk superconductor at B = 0. For superconducting films with 
thickness d much smaller than coherence length the role of the parallel component of the magnetic field is negligibly 
small. Thus, we can take account only of the magnetic field component B z perpendicular to the film surface and 
also neglect the dependence of the order parameter on z. For the sake of simplicity we restrict ourselves to the 
consideration of the one-dimensional case: B z (x) = H + b(x), where if is a uniform external magnetic field and b(x) 
is the z-component of the field induced by the magnetization M = M(x)zq (see Fig. 1). 
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FIG. 1. Superconductor/ferromagnetic (S/F) bilayer. 

Choosing the gauge A = A(x)yo, one can easily see that the momentum along the y axis is conserved, hence we 
can find the solution of the Schrodinger-like equation (1) in the form ^(r) = fk(%) exp(— iky), where function fk(x) 
should be determined from a solution of the one-dimensional problem: 

d 2 f k , {2TT V 1 

■jjT-Aix) -k) f k = -pTTT^fk ■ (2) 



dx 2 ' \<S>o" y ~ J "J JK £ 2 (T) 

Nontrivial solutions of Eq.(2) exist only for a discrete set of temperatures T n (k). The superconducting critical 
temperature T c should be define as the highest value max{T„(fc)}, corresponding to the lowest "energy level" 1/£ 2 (T) 
of the Schrodinger-like equation (2). 



II. SUPERCONDUCTIVITY NUCLEATION AT A DOMAIN WALL: AN ISOLATED ORDER 

PARAMETER NUCLEUS 

Let us start from consideration of a superconducting nucleus at a single domain wall taking the magnetization 
M near the wall in the form: M = Msign(x)zo (we assume that the domain wall width is much less than the 
superconducting coherence length). 



A. Domain wall in a thick ferromagnetic film: step-like magnetic field profile 

As it was mentioned above, for a rather thick ferromagnetic film (D 3> w) the expression for the distribution 
of magnetic field near the surface reads: B z = 47rAfsign(a;) + H, where H is an external applied magnetic field. 



2 



We choose the gauge in the form: A = (47rM|a;| + Hx)yo- At high temperatures the superconductivity far from 
the domain wall can be completely suppressed due to the orbital effect. On the contrary, near the boundary the 
superconducting nucleus can be still energetically favorable due to the mechanism analogous to the one responsible 
for the existence of H c3 critical field for superconducting nucleus near the superconductor-insulator interface (see, e.g., 
Ref. 24 ). Thus, a change of the magnetization direction which occurs at a domain boundary is responsible for a partial 
decrease of the orbital effect which provides conditions for the formation of localized superconducting nuclei at the 
domain walls at high temperatures (above the critical temperature far from the walls). Such a localized nucleus can 
appear only if we take account of proximity effect, i.e. consider the Cooper pairs to exist on both sides of the domain 
boundary. Such systems can reveal an interesting behavior in an external magnetic field. An external magnetic 
field applied to the sample results in a partial compensation of the field above one of the domains. As a result, 
the critical temperature of superconductor can depend nonmonotoneously on the applied magnetic field. Both the 
critical temperature of superconductivity nucleation far from the domain wall and critical temperature of formation of 
localized superconductivity at the wall should increase up to the external field value equal to the magnetic induction 
induced by the ferromagnetic moment. 

It is convenient to rewrite equation (2) in the following dimensionless form: 

-^ + (\t\ + ht-t ) 2 f k = Ef k , (3) 

where t = x/L, t = kL, 1? = $ /(2irB Q ), h = H/B , E = (T c0 -T) / AT° rb , the value AT° rb = T c0 $/L 2 characterizes 
the shift of critical temperature due to the orbital mechanism, and Bo is the maximum absolute value of the field b 
(in this subsection B = inM). 

For the case \t \ ~ * oo a superconducting nucleus will appear far from the domain boundary at a certain T c °°. In 
this limit the lowest eigenvalue E = |1 — \h\\ of equation (3) and, hence, the critical temperature is not disturbed 
by the presence of the domain boundary. On the contrary, for finite t values the superconducting nuclei to the 
left and to the right from the domain wall can not be considered separately due to the proximity effect. Provided 
the lowest energy level in the resulting potential well in equation (3) is minimal for a certain finite to coordinate, 
we get a superconducting nucleus localized at the domain boundary for temperatures above T c °°. The mechanism 
resulting in the appearance of such localized nucleus is analogous to the one responsible for existence of the surface 
superconductivity at the superconductor/insulator boundary for magnetic fields H c2 < H < H c3 . Indeed, for h = 
the potential well V(t) in Schrodingcr equation (3) is symmetric (V(t) = V{— t)) and the eigenvalue problem (3) can 
be considered only for t > with the boundary condition f' k (t = 0) = 0. For this particular case the energy minimum 
corresponds to tf t = E m i n = 0.59010 (Ref. 24 ). An increase in the h value will obviously result in an increasing 
asymmetry of the well V(t), and, thus, in the suppression of superconductivity localized at the domain wall. The 
equation (3) can be solved exactly in terms of Weber functions (see Ref. 24,28 ): 



A _< w (,/i+ft- 7 £ TI f s ), f >0 



(4) 



(5) 



Here C\ and C2 are constants, and the Weber function W(s,e) is the solution of the following equation: 

-^- + s2w = ? W (6) 
with the boundary condition W(s — ► +00, e) — > 0. Matching these solutions at t = we obtain: 

a/1 + hW' ( ^= —\ a/1 - hW' ( 42= 



w 



( *2 E_\ w ( to E_\ 



This equation can be solved numerically which allows us to obtain the function E(to, h). The resulting dependence 
of the critical temperature of superconductivity nucleation on parameter h is shown in Fig. 2. 
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FIG. 2. The temperature dependence of the critical magnetic field for a S/F system with a thick ferromagnetic layer. The 
solid (dashed) line corresponds to the superconductivity nucleation at the domain boundary (far from the domain boundary) . 

One can see that the external field suppresses the localized superconducting nuclei and the superconductivity 
localized at the domain wall exists only at a relatively weak applied field. As we increase an external field the 
superconducting nucleus shifts away from the domain wall towards the region where the absolute value of the total 
magnetic field is minimal. For < \h\ < 1 the curve E(h) calculated numerically can be fitted by the following simple 
expression: 



In this subsection we proceed with consideration of another limiting case D <C w and consider the problem of 
superconductivity nucleation in the field of an isolated domain wall in a thin ferromagnetic film: B z (x,z = 0) = 
4Mt&n~ 1 (D/x) + H. Obviously, for rather weak external magnetic fields H < Bo (in this subsection the maximum 
value of the domain wall field at z — is given by the expression Bq = 2irM) the superconducting order parameter 
nucleates in the region near the point x where B z (xo) = 0. Provided the localization length I of the superconducting 
nucleus is much smaller than the characteristic length scale of magnetic field distribution, we can expand vector 
potential as 




(8) 



B. Domain wall in a thin ferromagnetic film 




Such a local approximation is valid if the following conditions are fulfilled: 



B' z (x ) 



£ < 1 



and 



I < x . 



(9) 



Introducing a new coordinate t = (x — x )/£ we obtain the dimensionless equation 




(10) 




(11) 



4 



e= i( 1_ £) 



(12) 



(13) 



The lowest eigenvalue of Eq.(lO) eo — 0.904 is achieved at Q ~ 0.437. For the critical temperature T c of supercon- 
ductivity nucleation we obtain: 



This expression is valid when 



T c p — T c _ eo 
AT° rb ~ 7T \2BqD 2 



sin 1/3 (iJ/4M) 



1/3 



sin 



4/3 



4gj 

2B 



cos(ff/4M) 



4ttML» 2 



5>n 



(14) 



(15) 



Note that close to T c q the upper critical field has an unusual temperature dependence: oc (T c o — T) 3 / 4 . 

As we increase an external magnetic field H the position of superconducting nucleus shifts from infinity to the 
domain wall at x = 0. For rather large fields H the nucleus appears to be localized at the domain wall. Thus, the 
behavior of the nucleus coordinate in an external field is an opposite to the one considered in the subsection II A. The 
critical temperature for high field H limit is given by the expression 



1 - \H\/B . 



(16) 



The simple asymptotical formulas given above are in a good agreement with our numerical simulations of Eq.(2) (see 
Fig.3). 




(T c -T c0 )/ATf 



FIG. 3. The temperature dependence of the upper critical field for a domain wall in a S/F system with a thin ferromagnetic 
layer for BqD 2 /$o = 25 (solid line). The dash line corresponds to the analytical expression (14) at low fields, the dotted line 
corresponds to the high field asymptotics (16). 

For numerical analysis of the localized states of Schrodinger-like equation (2) with an external magnetic field 
we approximated it on a equidistant grid and obtained the eigenfunctions fk{x) and eigenvalues 1/£ 2 (T) by the 
diagonalization method of tridiagonal difference scheme. The typical behavior of the ground state wavefunction is 
shown in Fig. 4. 
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FIG. 4. The typical behavior of the ground state wavefunction for a domain wall in a S/F system with a thin ferromagnetic 
layer (solid line). The magnetic field profile is shown by the dash line. The parameters are BqD 2 /$o = 25 and H/B = 0.24. 



III. NUCLEATION OF SUPERCONDUCTIVITY FOR A PERIODIC DOMAIN STRUCTURE 

In this section we consider the effect of interaction of Cooper pair wavefunctions nucleated at different domain 
walls. Surely such interaction is important only for temperatures close to T c0 > w), otherwize for a rather 

large domain size w >• £(T) the overlapping of superconducting nuclei above different domain walls is exponentially 
small. For the sake of simplicity we consider here the case w <§C D and take the step-like distribution of magnetic 
field induced by the domain structure with the period a = 2w: b(x) = £?osign(a;), for \x\ < w and b(x + na) = b(x), 
where n is an integer. The corresponding vector potential can be chosen in the form: A(x) = Bq\x\ for |a;| < w and 
A(x + na) = A(x). 

In the absence of external field the general solution of Eq.(2) meets the Bloch theorem: 

f kg (x + a) = f kq (x)e^ a , (17) 

where q is a quasimomcntum. The nodcless wavefunction of the ground state corresponds to the value q = and is 
an even function of x, and, thus, we obtain /j(.(0) = /^(w) = 0. So we conclude that the solution at zero external 
field is identical to the one describing the superconductivity nucleation in a superconducting film of the thickness 
10 in the uniform magnetic field Bo. Following Ref. 24 we can obtain the ground state wavefunctions and energy 
E = (T c0 — T c ) /AT° rb as a function of the momentum k. The behavior of the resulting dependence of E{k) strongly 
depends on the parameter w/L. Two different regimes could be realized: for small values of w/L < 2.5 there is only one 
minimum of E(k) at k = w/(2L 2 ), which corresponds to the superconductivity nucleation above the domain center. 
For larger values of w/L one obtains two minima with equal energies E at fc™" and k™ n (fc™ 11 + k™ m = w/L 2 ). For 
w/L 1 the coordinates of these minima k™ ln and fc™ m and minimum energy E approach the values corresponding 
to the ones for isolated domain walls (see subsection II A) . Depending on the fc-momentum value the superconducting 
nuclei appear either above the walls at x — na (for fc™ m ) or at x = w + na (for k™ ln ). Thus, the nuclei at neighboring 
domain walls do not interact within the linearized GL theory. The dependence of the critical temperature on the field 
Bo m a periodic domain structure is described by the formula: 

1 - T c/ T c0 ^E 

where the function F(z) coincides with that for a superconducting film in the uniform magnetic field Bq which is 
plotted, e.g., in Ref. 24 . For a finite domain thickness w the critical temperature T c {Bq) appears to be larger than the 
one for a single domain wall. This difference in T c becomes rather large for small values z ~ when the nucleus 

is not localized near the domain boundary. For large z values one can obtain: F(z) — ► z/1.69, which corresponds to 
the dependence T C (B ) for a nucleus at a single domain wall. 
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If we apply an external magnetic field H the Bloch theorem is no more valid and the solution fk(x) appears to be 
localized. The energy level E(k) becomes a periodic function of the momentum k: E(k + AttHw/Qq) = E{k). The 
behavior of the upper critical field and structure of superconducting nuclei are controlled by the parameter w/L. The 
results of our calculations carried out using the same numerical scheme as in subsection II B are shown in Fig. 5. 




FIG. 5. The temperature dependence of the upper critical field for a periodic domain structure in a S/F system with a thick 
ferromagnetic layer for tyBqw 2 /$o = 5 (solid line) and nBow 2 /&o = 1 (dash line). 

For large values w/L the phase transition line is very close to the one found in the subsection II A, except for the 
small temperature region close to T c0 : AT <~ 4T c0 £o / w 2 ■ Outside this narrow temperature interval (and for H < B ) 
the wavefunction is localized at the domain walls (see Fig. 6). 
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FIG. 6. The behavior of the ground state wavefunction (solid line) localized at a domain wall in a periodic domain system. 
The magnetic field profile is shown by the dash line. The parameters are 7r_Bow 2 /$o = 5 and H/Bq — 0.095. 

The coordinates of these localized nuclei shift at ma as we change the momentum at 4:irHwm/<& { m is an integer). 
Let us note, that for rather weak magnetic fields H < B we observed a very peculiar behavior of the order parameter 
for a discrete set of field values given by the condition fc™" — fc™ m = 4irHwm/$o: the ground state wavefunction 
fk{x) has a two peak structure (see Fig. 7). 
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FIG. 7. The two peak structure of the ground state wavefunction (solid line) for a periodic domain system. The magnetic 
field profile is shown by the dash line. The parameters are ivBoW 2 /&o = 5 and H/B — 0.16. 

This fact is a natural consequence of the equivalence of the momenta k and k' = k + AnHwm/^o and the resulting 
resonant interaction of nuclei localized at domain walls separated by the distance w(2m — 1). 

For not very large values w/L < 2.0 the critical temperature becomes a monotonic function of the external mag- 
netic field because of the strong overlapping of wavefunctions corresponding to different domains. Therefore, the 
wavefunction is no more localized in a single domain (see Fig. 8). 
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FIG. 8. The behavior of the ground state wavefunction in a periodic domain structure for tyBow 2 /$o = 1 and H/Bo = 0.025 
(solid line). The magnetic field profile is shown by the dash line. 

However, even in this case we still observe a change in the slope of the phase transition line (see dash line in Fig. 5). 

The behavior of the upper critical field discussed above is not specific for step- like field distributions. To demonstrate 
this fact we studied the superconductivity nuclcation for the field profile B z {x) = Bo cos(2iTx/a) + H. The phase 
diagram on the plane H — T appears to be qualitatively similar to the one shown in Fig. 5. The critical temperature 
is a monotonic function of the external magnetic field for a/L < 4.5. For large parameters a/L » 1 and H < B 
the behavior of the critical temperature can be analyzed analytically following the approach used in subsection II B. 
The characteristic size of a superconducting nucleus and the critical temperature of superconductivity nucleation are 
given by the expressions: 



8 



T c \ Co 2 3 2n*B a* ( H 2 



(18) 



The validity range of this approximate description is defined by the conditions: 

H/B„ Y'^(2fB^V' 3 and f^X" » , . 



IV. CONCLUSIONS 



To summarize, we investigated the conditions of nucleation of localized superconductivity at the domain boundaries 
in hybrid S/F systems. The appearance of these localized superconducting nuclei should result in a broadening of 
the superconducting transition probed, e.g., by the resistivity measurements. We predict different regimes for the 
temperature dependence of the upper critical field near T c . The crossover between these regimes could be easily seen 
on experiment. In fact, the beginning of the resistivity decrease with the temperature decrease would correspond 
to the domain wall superconductivity, while its complete disappearance would signal the bulk superconductivity. 
External magnetic field would shrink the region of the domain wall superconductivity. Let us discuss some estimates 
of the physical parameters for the systems where the nucleation of superconductivity at domain boundaries could be 
observable. We can take, for example, the parameters of Nb (T c ~ 9 K and dH C 2/dT ~ 0.5 kOe/ K) and typical values 
of magnetization for ferromagnetic insulators 4nM ~ 1 — 10 kOe. The resulting increase in the critical temperature 
above a domain wall is quite strong: 8T C <~ 1 — 3 K. The thickness of a superconducting film must be much smaller 
than the distance between domains and ideal conditions correspond to the thickness of the order of several coherence 
lengths. So we conclude that the effects discussed above may be easily observed and could be quite important. Note 
that the behavior observed in Ref. 21 for S/F bilayers with bubble domains in a ferromagnetic film is qualitatively 
similar to our predictions. Generally the temperature behaviour of the critical field in S/F structures can be very 
rich (see Figs. 2, 3, 5) and it is strongly dependent on the domain structure and method of determination of the critical 
field. Careful measurements of the resistive and magnetic transition (including the measurements of the transition 
broadening) on the samples with a controllable domain structure would be very useful for the interpretation of the 
phase diagram and could give an important information on the domain wall superconductivity. 

Note in conclusion, that the existence of localized superconducting channels near the domain walls in S/F het- 
erostructures can provide an interesting possibility to realize a switching behavior provided we can move the ferro- 
magnetic domain wall. The superconducting channel in this case should follow the motion of the domain wall, which 
provides a possibility to control the conductance between certain static leads. 
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